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The bimetric variational principle is a subtle reinterpretation of general relativity that assumes 
the spacetime connection to be generated by an independent metric. Unlike the so called Palatini 
formalism that promotes the connection into a fundamental field, the new variational principle 
' results in a physically distinct theory since the potential for the connection carries new degrees of 

^ I freedom. The connection-generating metric naturally allows also an antisymmetric component. This 

sets torsion propagating! It is also shown here that while in the most straightforward generalization 
fT^ ■ of the Einstein-Hilbert action the nonmetric degrees of freedom become ghosts, there exist very 

simple actions which give rise to viable theories at the linearised level when subjected to the bimetric 
variational principle. However, the non linear interactions might bring unpleasant features like the 
Boulware-Deser ghost. This remains to be explored since this new type of bimetric theories does 
not, in principle, lie in the class of usual bimetric theories where non-linear interactions inevitably 
come in with new ghost-like degrees of freedom. 
U . 

PACS numbers: 04.20 Fy, 04.50 Kd, 98.80 Cq, 02.40 Ky 
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^ ; I. INTRODUCTION 

o| 

, Gravity is unique among the fundamental interactions of Nature because, while electroweak and strong interactions 
occur between particles in a background spacetime, gravity describes the dynamics of the spacetime itself. The 
] relation of gravity with the geometrical properties of spacetime lies at the heart of one of the most successful theories 
04 ■ of physics, General Relativity (GR). An issue of paramount importance is then to identify the degrees of freedom 
' ^ I [ associated with the spacetime geometry. 

J> . In GR, the dynamics of spacetime is described solely by its metric. However, from differential geometry we know 
; ""j ■ that spacetime, when described as a manifold, can be endowed also with an affine structure. To formulate physical 
rN . theories of gravitation, it is indeed necessary to introduce both metric and affine structures [l[ . Firstly, the need for 
■ measuring distances between two spacetime points^ or, equivalently, for measuring angles between vectors, forces us 
to introduce a metric structure. The affine connection is then needed to parallel transport vectors from one spacetime 
point to another. The concept of parallelicity is in principle independent of the metric structure. In GR however, the 
following relation is postulated: 

r"^j^ = ^g"^ {9p\,-i + 9\i,t3 - 5/37,a) ■ (1) 

Thus the connection is set to be the Levi-Civita, which clearly is a very special choice. It is symmetric and metric- 
compatible, i.e. Vaj^i/ = 0. When the latter condition is violated, we have non-metricity, which in general corresponds 
to D^{D + l)/2 independent degrees of freedom in a ZJ-dimensional spacetime. If the connection has an antisymmetric 
component, we are said to have torsion, which adds D^{D — l)/2 independent degrees of freedom. 



'Electronic address; |Jose.Beltran@unige.ch; agolovnev@yandex.ru; mindaugas@ugr.es; t.s.koivisto@astro.uio.no| 

^ Actually, the existence of a metric in a space is a stronger condition than the existence of a distance function, but wc shall assume that 
such distance function is the one compatible with a scalar product. 
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Recently it has been proposed that the spacetime connection could be generated by an independent metric 0]. 
Then, if we call this new metric g^j^, it acts as the potential for the spacetime connection as 

r V = ^5"^ {gpxn + 9>^i.P - 9M.>~) ■ (2) 

In the C-thcories of gravitation 0,01 8- conformal relation is assumed between the two metrics, g^j/ = C(72,)g^i/, where 
TZ ~ g^^'^R^viV) is the curvature of the spacetime. The C-theories introduce completely new theories and unify Einstein 
gravity and the so called Palatini gravity which correspond to specific choices of the conformal relation C (TV) . To sum 
up, Einstein's GR imposes the a priori rule ^ for the connection, and the C-theorics arc a conformal generalization of 
this rule. An alternative approach is not to prescribe a rule for the connection, but consider the connection-generator 
(jfj^i, as an independent tensor field. This approach we call the (unconstrained) bimetric variational principle. There the 
conformal degree of freedom adds effectively a massless scalar field 0] . This already shows that the resulting theory 
is physically distinct from metric GR, but what happens to the other degrees of freedom remains to be clarified. This 
is what we pursue in this study. 

However, it is first useful to briefly comment upon the difference of our new approach to the standard mctric-affine 
variational principle Q. The latter, often called the Palatini variation, promotes the connection into a fundamental 
field. It is well known that for the Einstein- Hilbert Lagrangians coupled to nonspinning matter, the "Palatini's 
device" results in GR^ , whereas for many classes of Lagrangians it generates underdetermined field equations 0, Q . 
More recent studies have discovered more problematical aspects of the resulting theories'^ [1343 , though they have 
some attractive properties in view of cosmology and even quantum gravity phenomenology, see the review [l^ . From 
different viewpoints, the culprit for the problems has been identified as the nontensorial nature of the connection field 
3 , the mixing of the metric-affine and Riemannian frameworks [l^, [l^l and the lack of dynamics • Adopting the 
bimetric variational principle removes all these possible issues. 

In section [n] we will derive the connection from a general tensor potential up to the leading order. In section Hill we 
shall analyze the antisymmetric sector of the theory. There appears an interesting model with propagating torsion. In 
section HVl we analyze the symmetric sector of the theory, and find that when applied to an action that is linear in the 
curvature, there appear additional ghost degrees of freedom. This rules out the bimetric variational principle when 
applied to the simple Einstein-Hilbert -like action. This leads us to consider a deformation of this action. Indeed, in 
section |V] we show that, by adding a standard Einstein-Hilbert term, the transverse traceless ghosts can be exorcised, 
although then the conformal mode becomes problematic. We conclude in section [VT] 

Throughout the paper a "mostly plus" signature of the metric is used (— , -|-, 



II. THE BIMETRIC VARIATIONAL PRINCIPLE 



Our starting point is to imagine such that we need not have Tj^j^j = 0, while according to ^ we have 



Vaff^^ = . 



(3) 



The field y^j, should be nondegenerate in order to have a well-defined inverse g'^". Since this field is nothing but the 
potential for the geometric connection, we do not need to impose further constraints on it. In particular, since it does 
not correspond to physical distances, it does not need to be symmetric 7^ .g,y^. We are thus led to consider the 
action 



S = 



-9 



1 



16nG 



y 



a 1 9^^\ 



(4) 



where ipa denotes the matter fields. Apart from a couple of hats added, this action is just the usual Einstein-Hilbert 
action. However, the differences are profound, as we shall find in the following. In order to proceed further we shall 
split the metric potential into its symmetric and antisymmetric components, namely: 



(5) 



^ To be precise, one needs to assume further conditions on the connection to uniquely recover GR in the Palatini formalism (sec for 
instance 0]). 

^ Some of the problems are cured when one considers a hybrid metric-Palatini theory Q . 
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with g^^ = 5(^1/) and S^j/ = g[fj,u] and we assume that they are not degenerate. Then, if we insert this decomposition 
into ([3]) we obtain the two sets of equations: 

^a9^.. = 0, (6) 
VaS^, = 0. (7) 

These equations determine the connection T^^ in terms of and E^^. Now, in order to give a more intuitive view 
of the theory we shaU also spht the connection as 

r^^ = r^^ + (8) 

with T^^ — r"^^) and S^^ = Tj^^j • Then, Eq. (jG]) and ([7]) will give us F and S in terms of g and S. This completes the 
exposition of the bimetric variational principle. It is clear we can decompose the fundamental degrees of freedom in 
gravity in this set-up to the usual metric plus the symmetric and the antisymmetric pieces of the connection-generating 
potential. We shall investigate each of the latter sectors separately. 

In the following we shall restrict to perturbative order. A perturbative solution to the equations above is given by 

f "7 = ^.9^" {d^9p\ + dfsgx^ - dxgp^) , (9) 
sip - i^^" (V^S.p + Vpl]^,„ + Vc«S/^p) . (10) 

where V is the covariant derivative for the connection T. This perturbative solution has been obtained at first order 
in S, assuming that it is small as compared to the symmetric part. 

The bimetric variational principle proposed here differs from other bimetric theories of gravity in several aspects. 
Concerning the theories in which one of the fields play the role of background metric and it is non-dynamical, the 
difference is clear, since we have that the two metrics are dynamical. Thus, our variational principle could only 
be equivalent to bimetric theories in which the two fields are dynamical. Moreover, one of our fields carry a non- 
symmetric component, so it could only be equivalent to bimetric theories with torsion. In any case, as we shall see 
later, when the antisymmetric part of g^u is small, it decouples from the symmetric component at quadratic level and 
one could expect the resulting action for the symmetric component to be equivalent to existing bimetric theories. To 
show that this is not necessary the case, let us consider the usual theories with two dynamical metrics /^jy and g^i^ 
that interact through a certain potential described by the action: 

S[.f^,u,g^u\ = ^ j d'^x^J~^R[f) + ^ j d'^x^/^R{g) + Sintif,,^, g,j.u]- (H) 

In this theory, the interaction term contains non-derivative couplings between the two metrics. However, in our pro- 
posed bimetric variational principle, the coupling between the two spin-2 fields is only through derivative interactions. 
In particular, it seems apparent that one cannot recover the Fierz-Pauli action for a massive spin-2 field from our 
variational principle for the linearized theory as one would require for Sint- We shall explicitly show this below. Thus, 
we do not expect the results obtained in those theories to be directly applicable to our case. 

Another difference appearing in our bimetric variational principle is that the field only appears as an auxiliary 
field or, in other words, it is a tensor potential for the connection of the spacetime. In that sense, this field is not 
physical and only the curvature generated by it can have physical effects, i.e., it can only appear through its associated 
curvature. This is crucially different from other bimetric theories because there both metrics are physical. This can 
be used as a guiding principle to add more terms in the action since g cannot appear directly, but only through R'^jj^s 
and its contractions with the spacetime metric g^i/. 

III. THE ANTISYMMETRIC SECTOR 

An interesting consequence of the restored dynamics is propagating torsion. At the classical level, torsion proves to 
be essential for total angular momentum conservation when intrinsic spin angular momentum is relevant (for reviews 
on torsion, see [l9l [20j). On the other hand, at the quantum level, a covariant version of Dirac equation also requires 
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the presence of a torsion. This necessity roots on the mass-reversal symmetry ^ of the Dirac equation when we want 
it to be a local symmetry, as required by special relativity. This could arguably be seen as one of the indications that 
torsion must be present in a fundamental theory of gravity, as discussed in [21| . As we will show below and contrary 
to the metric-afhne variation^, the bimetric variation procedure naturally yields a dynamical torsion field. 

Let us now consider the decomposition in ([8]). As one might expect, the symmetric component of the connection is 
nothing but the Levi-Civita connection of 5^,^, i.e.,® ^ag^i' = and the torsion field is determined by the antisymmetric 



component of the metric potential. An important constraint satisfied by S] 



Hi 



is that it is traceless 5*^^ 



0, as can be 



directly seen from the above expression. Then, it is straightforward to write down the Ricci tensor in terms of 
and 5?: 



5 A qol 



(12) 



where we have used the traceless property of the torsion tensor. Thus, when we plug this expression into the action, 
we obtain the following: 



S 



1 



167rG 



(13) 



Once again, this action might seem not to be new in the sense that it reproduces the usual Einstcin-Hilbcrt action in 
the presence of torsion. However, there is a crucial difference. The torsion field itself is not the fundamental field, but 
it must be expressed in terms of g and E, since it is them where the antisymmetric part of the spacetime connection 
emerges from. Thus, the fundamental fields of this theory are 5^1/, and 'Sap- At first sight, the true spacetime 
metric tensor t/^i/ might seem not to be dynamical because there is no kinetic term for it. This, however, turns out to 
be wrong because the equations of motion for g^^ will actually give evolution equations for g^^. Another interesting 
feature of the above action is that, since S^^ depends on first derivatives of Sa^, we indeed obtain a propagating 
torsion field even for the simple Einstcin-Hilbcrt term. Indeed, Eq. pop can be recast in the following form: 



7/3 



(14) 



where Fap^y is the strength field tensor of Eq,^ (Kalb-Ramond field), so that the gravitational action can also be 
written as: 



S 



d'^x. 



g' 



16nG 



1 



(15) 



where the dynamical origin of the torsion field becomes apparent. If we assume that the antisymmetric component 
of the connection potential Eq,^ can be derived from a potential, i.e., E = dcr, then F = and we obtain a gauge- 
fixing-like term for E. However, if we promote a into the fundamental field instead of E, we will obtain higher 
order equations of motion that could potentially introduce Ostrogradsky instabilities. Another possibility would be 
to make the torsion field be proportional to the Kalb-Ramond field so that the effective action becomes nothing but 
the standard Faraday action for a 2-form. The stability of the antisymmetric part may indeed require us to impose 
this constraint [1^. This could for instance be achieved by imposing a gauge symmetry on the torsion sector, as it 
happens for instance when one assumes that the torsion is fully derived from a tensor potential 
stability of each particular possibility should be carefully studied. 



24|. In any case. 



IV. THE SYMMETRIC SECTOR 



Next we study the stability of the symmetric piece. For this purpose, in the following we shall consider the perturbed 
theory around GR so that: 

gtj.1^ = gfj.v + h^iy (16) 



^ The mass reversal symmetry of the Dirac equation for the spinor field with mass m consists of the combined transformations: 

ip — > e'"'' 1/) and m — > — m. 
^ The metric approach can be rendered equivalent by a slight change of the Lagrangian [23 | . 
® The reader should be careful not to get confused between V and V. 
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with h^ii, a small perturbation. Since we focus on the symmetric sector in this section, let us neglect the contribution 
from S^^. This allows us to extend some of the results beyond linear perturbations. 
To proceed, we first notice that the inverse metric can be found as 



[g + hy 



n=0 



(17) 



where (—17 ^/i) = I, the unit matrix, and demanding that the series be convergent gives a precise meaning to the 
smallness of h. 

Now, the connection coefficients read 



'^^ap + .9'"' {daKp + dpKa " d^h^p 



(18) 



At the 0-th order in h it gives, of course, T^^. At the n-th order, n > 1, one has to use the n-th term of the sum with 
the first term in the brackets on the right, and the (n — l)-th term in the sum with the second term in the brackets. 
It yields 



^ a/3 — ^ aP 



r'^ 4- - 



^0 

E (-.9-^^) 



ig'"' {daKp + dph^a - dyhap) - g'^'^K-y^'lp 

gP" {VaKp + VpKa - '^Jlap) ■ 



An interesting feature is that 



1 



-.P - 2 

where the matrix g should be regarded as a function of matrices g and h. From the last expression we see that 

1 
2 



(19) 
(20) 

(21) 
(22) 



due to the symmetry of the metric g'^'^ . It also holds true at the linear level in h because the zeroth order approximation 
to g is equal to g which is also symmetric. Note that, generically, if we approximate the series for g only by terms 
with ^ n ^ A^, then the sum would not be symmetric. However, its asymmetry is always next order in h. Indeed, 



N 



E (-^-^^)' 



n=0 



E i~9~'h) 

n=0 

I-{-g h) 



E i-a-^'^y 

.n=N+l 



.n=0 



Next, it is necessary to evaluate the Ricci tensor, 

p pi pa pi pa I pa y-p pa y-p 



(23) 
(24) 

(25) 



Substituting r"j^ = + dV^^, where the tensor quantity SV^^ is already known from eq. (jl9p in the form of 
convergent (for small enough h) power series, we easily get 



(26) 
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The derivative part of the Ricci tensor variation is purely a surface term, and can be neglected. ^ Then, the resulting 
action is: 



S = 



1 



16nG 
1 



1 



IGttG 



n=0 



Kp 



16nG 

1 /3 



d'^X\/—gR 



(27) 



We see that the added ingredient is just a symmetric tensor field h possessing a quadratic in V/i action, although 
with non-linear in h non-derivative coefficients. At the first (quadratic) order the coefficient in front of the covariant 
derivatives reduces simply to jg^^g^"^ ■ (In general, this coefficient is equal to jg°'''g^'^ as a function of g and h. And 
this is actually a non-perturbative result as long as E^j, vanishes, valid at any, not necessarily small, values of h.) In 
this limit the action in eq. ([27|) is just equivalent to the linearized Einstein-Hilbcrt action: 



S 



1 



IGttG 



(28) 



where we have assumed that the metric g^^ is the Minkowski one, and therefore all covariant derivatives are substituted 
by the ordinary ones. 

One of the two gravitons has inevitably ghost excitations. Indeed, we were assuming the diag(— , +, +, +)-signature 
of the metric throughout the paper which implies that the /i-field is a tensor ghost. Of course, if we switch to the 
opposite sign convention, then the /i-field is healthy. But after separating it, we have the GR action oc J d'^XyJ—gR 
while, with the new convention, it must have an opposite sign. Let us demonstrate the problem explicitly by expanding 
the metric around the Minkowski spacetime g^u = rjp_i, -\- hp^^. The second order action for h acquires the form 



1 



167rG 



d^x 



(29) 



The first term we already know from eq. ([^5]). S'^'^'^Rpi, ~ '^r^^'^r^j, — ^^'^pST^^ -f surface terms, and it gives precisely 
the same action for h as the action for h in eq. ([28l) . However, it is not yet the end of the story. With S'^^^Rpi, ~ 
daSVil^ - ^^'^r;^^ and 5'-^'> {^^g'"') = -h'"' + i/i^??'"'', we write the second term in 1^ in the form 



fict u 



which, after some integration by parts, transforms to 
and, together with the first part, it gives the final result: 



5(2) 



1 



647rG 



dc^hp^d'^h^" - 2d''hp,d''h^^ + 2dc,h''^'dph$ ~ dphy^h 



(30) 



(31) 



(32) 



^ Nevertheless, those derivatives can also be presented in a nice convergent form. The only ingredient we need is the covariant derivative 
of the power series 

oo oc 

where the f)-matrix simply denotes the hJJ-field. It is easy to see that differentiating the (— fl)" term one gets n terms with — Vfl in all 
possible positions. If we differentiate the whole sum, but consider only the terms with (— f))™ to the right from — V(), then there would 
be exactly one possible term with any fixed power of — () to the left, where (—1))* comes from the differentiation of (— (jJ'^+'^+l term in 
the initial sum. Therefore, it is easy to see that 



(V^/i?) 



EM)* 



Obviously, this formula generalises the differentiation rule for a function f{x) = (1 + g(x)) ^ . 
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We see that the tensor modes are ghosts, where TT stands for spatially transverse and traceless. 

The total second order action for both types of perturbations acquires the form 



5(2) 



1 



167rG 



(33) 



At this level we have two separate diffeomorphism invariances (one of which is fake), and therefore can not judge 
upon the scalar and vector degrees of freedom. However, we definitely see that the two types of tensor modes do have 
opposite signs of kinetic functions. 

In summary, the theory is full of ghosts, at least in the simplest version of this class of models which we have just 
considered. After all, it is not probably that much surprising. Let's forget about the tensor nature of the metric fields. 
We take two scalar variables instead, x for g and y for g. Suppose we have an action of the form S = J dt{x + x^)y, 
somewhat analogous to our gravity model. Integrating by parts we have S = J dt{—iy + i^y). After diagonalising 



4^ 



where z = logy, and see that one of 



this action we get S" = J dt y [{x - \ - \ = J dt 
the modes is always a ghost. Finally, it is interesting to note that no mass terms appear at this level for none of the 
metrics, as it was expected from our discussion above pointing out the differences of this theory with usual bimetric 
theories. The reason for the absence of the mass term is that, even though the fields will interact at higher orders, all 
the interactions will be through derivative couplings. 



V. ON MORE ELABORATE THEORIES 



To render the theory viable, we have to reconsider our action. A very simple way to extend the previous version of 
the theory is to consider the curvature TZ = g^^^R^i, to be added to the usual Einstein-Hilbert term R = This 
way it loses some of the elegance, but wc can use it as a proof of concept. Let us say the relative coupling strength of 
the nonmetric degrees of freedom is given by the parameter < k < 1. This adds one parameter to the theory that 
can then in principle be constrained by the experimental tests of gravity. So, we write the action as 



S = J d^xy^g^'' (R^, - kR^^j . (34) 

Using the calculations of the previous section it is easy to see that, at quadratic order and with Newton's constant 
renormalised by G — > G/ (1 — k), the action becomes: 

5(2^ - --^ I d^x[(d^h^,d^h^^'' ^2d"h^^d''h^^ + 2d^h^''dyp-d^hy^h^p^ 



L hi \ / ± Ki 

where now the two graviton fields arc healthy as long as^ < k < 1. That is, wc obtain two copies of a massless 
spin-two fields plus the 2-form field determining the torsion of the geometric connection. It is remarkable that 
perturbations in both symmetric fields h^^ and h^^i, happen to acquire the usual kinetic term of linearised GR, which 
is well-known to be precisely the only viable action for pure massless spin-two fields^ [23j . Thus, at the leading order, 
wc have two sets of spin-2 fields entering into the action in a symmetric way. When matter is present, this symmetry 
is broken. Standard bosonic matter fields couple only to the spacetime geometry given by g^jy. This, for instance, is 
the case for scalar fields with standard kinetic terms or gauge (both abelian and non-abelian) fields constructed out 
of their strength fields"'^''. Of course, when non-minimal couplings are present, matter fields will also couple to g^,^. 
This will be, for instance, the case of scalar fields [2^ with higher order derivative terms, like K-essence or galileon 
fields, or vector fields with direct couplings to curvature. On the other hand, fcrmionic matter like spinor fields are 



* This is true for the signature convention with diag( — , +,+,+)■ 

^ Though recently the result was generalized to nonlocal theories [25{ . 

" This statement depends on how we decide to covariantize the definition of the strength field. Let us consider the electromagnetic field 
defined in flat spacetime as F/n, = dfj,Ai, — d^Ai^. If we just replace partial derivatives by covariant derivatives, _F^y will acquire a term 
involving torsion. However, if we define the strength field as the exterior derivative of we do not have couplings to the connection. 
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nonminimally coupled to the spacetime connection, and thus wiU always act as a source also for g^u. Since a Dirac 
field only couples to the totally antisymmetric part of the torsion field, such a coupling will be mediated by F^p^, 
which is defined after ([T4|) . Of course, fermionic fields with higher spin will couple to the other components of the 
torsion field. Notice that, whereas nonminimal couplings of bosonic fields will lead to classical effecs, fermionic matter 
can only give rise to effects at the quantum level. However, since the geometric connection has a non-symmetric 
part, i.e., torsion, one could couple this connection to fermionic fields, whereas bosonic fields would only couple to 
the matter connection. Moreover, the considered distinction between the metric and affine structures of the physical 
spacetime resolves the covariantization of the Dirac equation that leads to the need of a connection with torsion, as 
explained above. One may even speculate on making contact with supersymmetry by a possible relation of g^^ with 
the supersymmetric partner of the graviton that couples to fermions j27l |28| . 

To end, let us mention one important issue. The problematic part of quantizing Einstein gravity, and a big obstacle 
on the way of deforming it classically, is the conformal degree of freedom. There is indeed a scalar mode with the 
negative sign of kinetic function. Fortunately for classical GR, it is killed by the gauge invariance. However, once you 
break it, you should expect some problems to come about, such as the sixth mode or the Boulware-Deser ghost. In this 
bimetric context, we encounter the mirror image of this problem. Since the conformal degree of freedom for the action 
(HI is a healthy scalar field [H, it seems that by curing the other nonmetric degrees of freedom in the theory (|34p . 
essentially by the minus sign in front of the coupling k, we have flipped the conformal mode into a ghost-like scalar. 
A simple way out would be a prescription where the conformal difference between the two metrics is not allowed to 
propagate^^. Whether this is actually necessary in order to construct theories that are completely ghost- free to any 
order in perturbation theory, is a problem we cannot address within the scope of the present study. 

VI. CONCLUSIONS 

The results presented here invite to study the nonlinear theory ensuing of bimetric variation. The main difficulty to 
go beyond the linear regime is to fully solve the equations ([3]) that give the connection in terms of the potential metric 
fields expressed by means of g^i, and S^j^ . Let us note that recently, in the context of ghost- free massi ve g ravity, the 
nonlinear couplings of two metrics that avoid the Boulware-Deser instability, have been discovered |29l - l33j . Whether 
we can obtain those couplings from the bimetric variational principle remains to be explored. Also, because of the 
natural presence of propagating torsion in the considered bimetric formalism, some relations with torsion massive 
gravity, in which the torsion field provides the graviton with a mass (s^. Issj. would be interesting to study further. 

To conclude, in this paper we intended to present a new variational principle in which the affine connection is 
derived from a general tensor potential, which is different from the spacetime metric. We have studied the linear 
regime and found that the theory is equivalent to GR plus an additional pure spin-2 field and a 2-form leading to 
propagating torsion. However, the healthiness of the additional spin-2 field forced us to consider a more complicate 
action (j34p than the one consisting merely of a term linear in the curvature TZ. Unfortunately, having cured the 
tensor modes, we rendered the conformal degree of freedom problematic. It remains to be seen whether it is possible 
to construct a completely viable theory in this formalism. In this respect, it is worth mentioning that the proposed 
bimetric theory is, in principle, different from those already existing in the literature so that the results about their 
stability cannot be directly applied to the present case. 
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